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ABSTRACT 

We simulate the evolution of a 10 9 M dark matter halo in a cosmological set- 
ting with an adaptive-mesh refinement code as an analogue to local low luminosity 
dwarf irregular and dwarf spheroidal galaxies. The primary goal of our study is to 
investigate the roles of reionization and supernova feedback in determining the star 
formation histories of low mass dwarf galaxies. We include a wide range of physical 
effects, including metal cooling, molecular hydrogen formation and cooling, photoion- 
ization and photodissociation from a metagalactic (but not local) background, a simple 
prescription for self-shielding, star formation, and a simple model for supernova driven 
energetic feedback. To better understand the impact of each physical effect, we carry 
out simulations excluding each major effect in turn. We find that reionization is pri- 
marily responsible for expelling most of the gas in our simulations, but that supernova 
feedback is required to disperse the dense, cold gas in the core of the halo. Moreover, 
we show that the timing of reionization can produce an order of magnitude difference 
in the final stellar mass of the system. For our full physics run with reionization at 
z = 9, we find a stellar mass of about 10 5 M at z = 0, and a mass-to-light ratio 
within the half-light radius of approximately 130 Mq/Lq, consistent with observed 
low-luminosity dwarfs. However, the resulting median stellar metallicity is 0.06 Z , 
considerably larger than observed systems. In addition, we find star formation is trun- 
cated between redshifts 4 and 7, at odds with the observed late time star formation in 
isolated dwarf systems but in agreement with Milky Way ultrafaint dwarf spheroidals. 
We investigate the efficacy of energetic feedback in our simple thermal-energy driven 
feedback scheme, and suggest that it may still suffer from excessive radiative losses, 
despite reaching stellar particle masses of about 100 M , and a comoving spatial 
resolution of 11 pc. 

Key words: galaxies: dwarf - Local Group - hydrodynamics - galaxies: star forma- 
tion - dark ages, reionization, first stars - ISM: evolution 



1 INTRODUCTION 

Dwarf galaxies are the most dark matter dominated ob- 
jects in the universe, with gravitational masses a hundred or 
even a thousand times greater than their observed baryonic 
masses. They therefore appear to lie at the extreme limit 
of galaxy formation. Despite this, however, they exhibit re- 
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markable consistency in many of their properties. Milky Way 
dwarfs appear to have a constant central mass (within 300 
pc) of 10 7 M over five orders of magnitude in luminos- 
ity (|Strigari et al. 20081) and their mean stellar metallicity 



tightly correlates with their luminosity over more than three 
orders of magnitude ( Kirby et al.|20lf ). Despite their small 
sizes and large mass-to-light ratios, dwarf galaxies appear to 
lie on many of the scaling relations seen in massive galaxies 
flKormendy et al.|2009| |Tolstoy, Hill, fe Tosi||2009| ) . Under- 
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standing the physics that shape dwarf galaxies is partially 
an effort to understand a particularly inefficient mode of 
galaxy formation, but it also is an effort to understand how 
physical processes involved in galaxy formation scale across 
galaxy mass to produce the observed continuous evolution 
in galaxy properties. 

The competition between a variety of baryonic heating 
and cooling mechanisms plays an important role in regulat- 
ing the gas available for star formation in galaxies. Under- 
standing galaxy formation is therefore in part about under- 
standing how these microphysical processes work in galaxies. 
These processes operate in a non-linearly growing potential 
well while being perturbed by environmental effects such 
as ram pressure and by effects of secular evolution of mas- 
sive baryonic structures such as disks and bars. Since dwarf 
galaxies in the Local Group, especially dwarf spheroidal 
galaxies (dSphs), are observed to have large mass-to- light 
ratios and old stellar populations, it is likely that their evo- 
lution has been heavily influenced by heating processes that 
suppressed later star formation. In addition, the number 
of dwarf galaxies observed around the Milky Way is sev- 
eral orders of magnitude less than the predicted number of 
dark matter subhaloes found in cosmological N-body simu- 
lations ( |Klypin et al.|1999| [Moore et al.|1999[ ). This 'missing 
satellite' discrepancy had been somewhat alleviated by the 



2000 



Joung & Mac Low 2006), but how these processes 



discovery of many 'ultrafaint' dSphs (e.g. Willman et al 



2005 ) and by hydrodynamical simulations of dwarf satellite 



systems that have shown baryonic processes can suppress 
star formation ( |Okamoto et al.||2010| |Wadepuhl &; Springel| 
2011). Several heating processes are at play in dwarf evolu- 
tion, including the metagalactic UV background during the 
epoch of reionizatio n (|Bullock, Kravtsov, &; Weinb erg 2000 ; 
Gnedin fe Kravtsov||2006), supernova feedback (Governato 
et al. 2010; Sawala et al. 2010) as well as local UV heating 
and winds from young stars. 

While these heating processes act in the evolution of 
galaxies of all sizes, dwarf galaxies provide an attractive lab- 
oratory for their study since their small sizes allow simula- 
tors the ability to reach higher resolutions. High resolution 
is important for resolving the dense gas from which stars 
form and for retaining the energy in hot gas, and is there- 
fore key to constraining the star-formation rate. High spatial 
resolution (of order 1 pc) is especially important in limiting 
diffusion from cold, high-density gas into hot, rarefied gas 



produced by supernova explosions (e.g, de Avillez &; Mac 
Low||2002| |joung fe Mac Low||2006| ). Without resolut ion on 



this scale, feedback energy is rapidly radiated (e.g. Katz, 
Hernquist, & Weinberg 1992). Simulations of galaxy forma- 



tion often compensate for this behavior by imposing feed- 
back driven winds or by enhancing the impact of feedback 
heating by either distributing energy to diffuse gas phases 



(e.g. Springel & Hernquist 2003 |Scannapieco et al. 2006) 



or turning off cooling for a time ( |Stinson et al.||2006 Gov 
ernato et al. 2010). These different approaches to modeling 
feedback can produce different star formation histories and 
galaxy properties even when keeping the total amount of 
supernova energy the same between different model imple- 



mentations (Schaye et al. 2010 Sales et al. 2010). We un- 



derstand a great deal about the structure of supernova rem- 
nants from observations of nearby resolved examples (e.g. 
Badenes 2010), and it is possible to model the impact of su- 
pernovae in idealized simulations of the ISM (e.g, de Avillez 



depend on, and interact with, global galaxy properties to 
produce global behavior such as galactic winds remains an 
open question. 

Heating mechanisms appear to be very important in 
dwarfs, but understanding the cooling mechanisms that bal- 
ance them is also crucial. Many of the heating processes 
we have outlined are a consequence of stellar evolution and 
are therefore correlated with the star formation rate. In low 
mass haloes with very low metallicities, the effect of molec- 
ular line cooling from H2 may be an important factor in 
determining the star formation rate. Cosmological galaxy 
simulations often assume an equilibrium cooling curve for 
low metallicity gas, but non-equilibrium effects may also be 
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important (Abel, Bryan, & Norman 2002 Glover Sz Clark 



As we have discussed, reproducing the proper phase 
distribution of gas in the interstellar medium is quite dif- 
ficult, however, modeling star formation from this medium 
through the creation of aggregate star particles is also quite 
uncertain. We see stars forming from dense molecular clouds 
in our own galaxy and Ke nnicutt] ([T998 ) has demonstrated 
that the star formation rate of galaxies correlates strongly 
with the gas surface density. Despite this observed rela- 
tion, the gas density threshold for star formation is often an 
adjustable parameter in cosmological galaxy simulations in 
part because of resolution limitations, and indeed, adjusting 
the star formation threshold at fixed resolution can produce 



very different galaxy properties (see Guedes et al. 2011). 
Modeling star formation from particular gas phases may also 



produce different behavior; |Kuhlen et al. ( 2012[ ) have shown 
that modeling star formation as dependent on the molecular 
gas fraction suppresses the star formation rate in low mass 
galaxies at high redshift. However, Glover & Clark ([2012) 
have argued that simply the gas density determines the star 
formation rate and the molecular gas fraction plays little 
direct role. Since many of the other uncertain subgrid pro- 
cesses that we have discussed (i.e. feedback) are dependent 
on the star formation rate, uncertainties in modeling star 
formation have the potential to compound themselves in a 
galaxy model, and in low mass dwarf galaxies where star for- 
mation and its quenching can balance so precariously, these 
issues may be especially important. 

For this study, we have deliberately chosen simple and 
widely used prescriptions to model subgrid baryonic effects 
such as supernova feedback and reionization. Our goal is to 
test these model prescriptions at high resolution by carefully 
comparing our simulation results to observations. We believe 
that the specific ways in which our simulations succeed or 
fail to reproduce observations can tell us much about the 
ability of these model prescriptions to capture the relevant 
baryonic physics properly. 

We will also attempt to better understand the factors 
that control the gas content, star formation rate, and stellar 
mass for Local Group dwarf analogues. We are particularly 
interested in probing the relative impacts of supernova feed- 
back and reionization on the evolution of these low-mass sys- 
tems. We will attempt to do this by better understanding the 
detailed evolution of the gas phase distribution within a cos- 
mological context. The simplicity of our model assumptions 
for subgrid baryonic physics will hopefully provide clarity in 
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interpreting the implications of our results for the evolution 
of real systems. 

To these ends, we have conducted a series of high reso- 
lution cosmological, hydrodynamical simulations of a single 
dwarf halo with prescriptions for H 2 and metal line cool- 
ing; star formation and supernova feedback; and a global 
UV background that ionizes neutral gas and dissociates 
molecules. We target a dark halo mass of 10 9 M . Mod- 
eling of the dynamics of Local Group dwarfs suggests they 
reside in dark matter haloes of approximately this mass or 
smaller (e.g |Tollerud et al.|201l) . 

We describe our simulation methods and physical pre- 
scriptions in Section [2] and present the results of our canoni- 
cal simulations in section |3J"1 We have also conducted simu- 
lations where we have adjusted some of the physical prescrip- 
tions implementing supernova feedback and reionization to 
gauge their effect in our model. We present the results of 
these alternate physics runs in Section [3^2] We discuss the 
results of our simulations and compare our results to obser- 
vations and other simulations in Section [4] and summarize 
our conclusions in Section |5] 



2 METHODS 

We conducted our simulations with the structured adaptive 



mesh refinement (AMR) hydrodynamics code Enzo (Bryan 
1999||Norman fe BryanllQQQllO^ea et al.|2004| ). Enzo uses 
an Eulerian method for solving the hydrodynamics equa- 
tions on a Cartesian grid. We chose a grid-based method 
because of its ability to resolve interfaces in multiphase gas 



( Slyz et al.|2005||Agertz et al.|2QQ7| |Tasker et aLp 008). We 
used the ZEUS hydro solver in Enzo which uses a simple 



and robust second-order finite difference method (van Leer 



1977 Stone & Norman 1992). Dark matter and stars are 



represented by particles, and the gravitational interactions 
are computed by solving the potential on the mesh ( |Q'Shea 
|et al.|[2004| . Our model also solves non-equilibrium evolu- 
tion rate equations for e~, H, H + , He, He + , He ++ , H~, H2 



and Hj, as well as cooling from these species ( Anninos et al 
1997 Abel et aL| [l997). In addition, we include metal line 
cooling, photoionizing and photodissociating backgrounds, 
star formation and thermal feedback, and self-shielding by 
neutral and molecular gas, which are all described in more 
detail below. Our version of the code is publicly available 
in the auto-inits branch of the online Enzo repository as 
changeset 15651fe320ff. 



2.1 Initial conditions 

We generated cosmological initial conditions with Enzo's ini- 
tial conditions generator inits. We used the best fit cosmo- 
logical parameters from the WMAP5 data release (Q m = 
0.274, Q b = 0.0456, Ov = 0.726 and h = 0.705) flHinshaw 



et al. 2009). We begin our simulations at redshift z = 99 



when perturbations in the universe still grow linearly. Our 
cosmological box is 4 comoving Mpc h _1 on a side and is 
divided into 128 3 cells on the root grid. We note that our 
box size is too small to accurately sample large scale modes, 
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and the largest mode in the box becomes non-linear at low 
redshift. This will affect the statistics of the haloes in our 
box, and will also impact accretion at late times; however, 
since this is an exploratory study to understand the gross 
evolutionary properties of a single low mass dwarf halo, we 
argue that this level of error is acceptable. In particular, we 
note that while the accretion history of our chosen halo may 
be biased, it is likely to fall within the wide range of possible 
such histories for haloes of this mass. We will redress this 
problem in later work. 

Our goal is to simulate a single dwarf halo at high res- 
olution. We first selected a target halo from a simulation 
of our cosmological box with no refinement. We identified 
haloes close to our target mass of 10 9 M at z = 0. We 
found that lower mass haloes traveled a longer distance rel- 
ative to their size than higher mass haloes (reflecting the 
fact that all haloes participate in large-scale bulk flows) . We 
chose a halo that traveled a relatively small distance rela- 
tive to its r2oo radius during the simulation to minimize the 
area of refinement necessary in resimulations. Our chosen 
halo travels a distance of 0.9 comoving Mpc h _1 during the 
simulation. This choice was motivated by concern for com- 
putational efficiency, however, it had the result of selecting 
a relatively isolated halo. The closest halo with comparable 
mass to our target halo has a mass of 8.52 x 10 8 M© and is at 
a distance of 338 kpc, which is approximately 14 times r2oo 
for our halo. The most massive halo in the box is 1.43 x 10 12 
M© and is 2.9 Mpc from our target halo. 

The target halo identified, we then created new, refined 
initial conditions based on the trajectory of our halo and 
its particles. We laid down a series of nested static subgrids 
that refine the root grid by a total of three additional levels 
so that in the most refined region the effective initial resolu- 
tion was 1024 3 . The cell size within the most refined static 
subgrid before adaptive refinement is therefore 5.54 comov- 
ing kpc. One corner of the maximally refined static subgrid 
was defined by the target halo's final position plus a buffer 
of 10 r2oo radii. The opposite corner was set to contain the 
volume enclosing the initial positions of the halo's particles 
in the coarse simulation, plus an additional buffer of 10% 
of the width of this enclosing volume in each direction. A 
buffer region of four intermediate resolution grid cells was 
placed around the boundaries of each successively nested 
static grid. The initial resolution of a cell sets the mass of 
the dark matter particles laid down within it, so the higher 
resolution cells have lower mass dark matter particles. The 
buffer regions are important to ensure that higher mass par- 
ticles did not impact the evolution inside the highly refined 
region. These initial conditions set the minimum mass of 
dark matter particles to be 5353 M©. 

To further restrict refinement to the halo of interest, 
we have employed a new technique that identifies dark mat- 
ter particles in our refined initial conditions that end the 
simulation in our target halo and then requires the pres- 
ence of these particles for the refinement of a given cell. We 
identified halo dark matter particles in our refined initial 
conditions by first running a dark matter only simulation 
to z = 0. We identified particles that ended the simulation 
within 3r2oo radii of the halo centre and tagged them as 
'must refine particles.' In all subsequent runs with hydrody- 
namics, we restrict adaptive refinement to cells that contain 
must refine particles and that are within the highest reso- 
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lution initial subgrid. Cells meeting these criteria may be 
further refined based on their dark matter and baryon den- 
sity such that additional refinement was added (by a factor 
of 2) whenever the dark matter mass in a cell exceeded four 
times the dark matter particle mass, with a similar crite- 
rion for the gas. We did not explicitly implement refinement 
based on the Jeans length, however, we found that with our 
density refinement criteria that in most dense regions we ei- 
ther resolved the Jeans length by four grid cells (meeting the 
criterion of |Truelove et al.' ([1997) to avoid artificial fragmen- 
tation) or the density of the gas was greater than our chosen 
threshold density for star formation. The net result of our 
method is that the adaptive refinement tracks the dark mat- 
ter particles and gas that ends up in the halo of interest. In 
our highest resolution run (Simulation RIO in Table [T}, we 
find that 99% of the dark matter particles within r2oo are 
must refine particles at z = 0. This method resulted in a 
computational savings of more than a factor of 10, as com- 
pared to allowing refinement anywhere in the most refined 
initial grid, which is standard practice. 

In our highest resolution simulations, we allow a total of 
12 levels of adaptive refinement beyond the root grid. These 
simulations therefore have a minimum cell length of 10.8 
comoving pc. This means at high redshift (z ~ 9) our calcu- 
lations have a physical resolution of order one pc. We have 
also conducted simulation R43 that allows only 10 levels of 
adaptive refinement beyond the root grid for a minimum cell 
length of 43.3 comoving pc. We will discuss this simulation 
in Section [4~3l 



2.3 UV background 

We include both photoionizing and photodissociating meta- 
galactic backgrounds in our model. We assume spatially uni- 
form photoionizing rates for H I, He I and He II that vary 
with redshift, taken from Haardt h Madau (2001). This 



background includes ionization from both quasars and stars. 
For most of our runs we start turning on the photoioniz- 
ing background at z — 7 and gradually ramp it up to full 
strength by z = 6. The photodissociation rate of H2 is taken 
to be 1.13 x 10 -8 times the flux in the Lyman- Werner bands 
( Abel et al.|1997 ). We adopt the average flux in the Lyman- 
Werner bands from the spectra of Haardt & Madau (2011 ). 



This flux is redshift dependent and we assume it to be spa- 
tially uniform. We turn on the photodissociating background 
at the same time as we start the photoionizing background 
(without a ramp-up period, as there is no significant gas 
heating associated with this transition). In order to investi- 
gate the impact of the timing of the reionization epoch, we 
also carried out a simulation in which we turned on the UV 
background at an earlier redshift, z = 8.9, and ramped up 
the ionizing background to full strength by redshift z = 8. 

Although the UV background is spatially uniform, 
we apply a crude self-shielding correction. The H I ion- 
izing flux felt by each cell is attenuated by a factor of 
exp(— AxnHKJion), where Ax is the cell size, nm is the neu- 
tral hydrogen density, and ai on is the frequency-averaged 
photoionization cross- section for H I. We computed the 
frequency- averaged photoionization cross-section as 



"0 

00 ' 



(i) 



2.2 Metal cooling 

As noted earlier, we include radiative cooling from nine 



atomic and molecular species ( |Abel et al.||1997 Anninos 



et al. 1997), whose (non-equilibrium) abundances are ex- 
plicitly followed. In addition, the code also tracks a mean 
metallicity, and additional metal heating and cooling rates 
are interpolated from equilibrium cooling tables, in a way 
similar to that described in |Smith, Sigurdsson, fe~ Abel 
(2008): Cloudy ( Ferland et al.|1998| is used to generate five- 
dimensional tables that give the net metal cooling and heat- 
ing rates, depending on the gas metallicity, electron frac- 
tion, density, temperature and redshift. The redshift tracks 
the time-dependence of the metagalactic UV background 
(see below). The combination of heating and cooling from 
non-equilibrium atomic and molecular hydrogen combined 
with equilibrium metal lines, allows for a more accurate es- 
timate of the cooling rates, particularly for low metallicity 
gas. Since we follow the abundance of H 2 , the cooling rates 
we compute for low temperature, low metallicity gas where 
H2 cooling dominates are particularly accurate. We do not 
include molecular hydrogen formation on dust grains, but 
by the time this pathway becomes important, metal cooling 



dominates over H2 cooling anyway (Glover 2003). Finally, 



we note that we do not include molecules formed out of 
heavy elements, but these tend to be less important than 
fine- structure lines ( | Glover &; C lark 2012), particularly at 
low abundances. 



where J is the flux of the ionization background, which we 
approximate as J oc ^, and a is the frequency dependent 
ionization cross-section, which is zero for frequencies less 
than 1/0, the threshold frequency for ionization of H I, and 
a oc A is~ 3 for frequencies greater than z/ . The parameter f3 
is the power law slope of the ionization background which we 
take to be -1.57, approximately the slope of the photoion- 
ization background that we use. The parameter A is the 
peak ionization cross-section for H I. With these assump- 
tions we calculate the average photoionization cross-section 
to be a m = A /3/(f3-3). 

A similar correction is made for the H I photoheat- 
ing rate, where we use the frequency- averaged photoheating 
cross-section which we compute as 



Cheat 



/ l£l{hv-hv )a&v 
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00 ' 

/ ^(fci/-fcm)di/ 



(2) 



which we find is proportional to ai on by a factor of (/3 + 
l)/(/3 — 2). Corrections are also made for He I and He II 
ionization and heating using the corresponding values of Aq 



for each species taken from Brown ( 1971 ) and Osterbrock & 
|Ferland ( 2006| . These corrections are approximate, as they 
assume that only gas in the local cell contributes to self- 
shielding, but they do suppress ionization and heating in 
dense regions. 
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Table 1. Summary of Simulations. 



Name 


(comoving pc) 






Az 


H 2 
Cooling 


(z 1 ©) 


RIO 


10.8 


3.7 x 10- 


6 


7-6 


yes 


none 


RIO-earlyUV 


10.8 


3.7 x 10- 


6 


9-8.9 


yes 


none 


R10-noH2 


10.8 


3.7 x 10- 


6 


7-6 


no 


none 


RIO-noUV 


10.8 


3.7 x 10- 


6 


none 


yes 


none 


RIO-noFB 


10.8 







7-6 


yes 


none 


RlO-noFB-LimCool 


10.8 







7-6 


yes 


0.1 


RIO-lowFB 


10.8 


io- 6 




7-6 


yes 


none 


R10-DM 


10.8 












R43 


43.3 


3.7 x 10- 


6 


7-6 


yes 


none 


Note: The quantities presented in each column 


are 


(1) the 


minimum 


cell 



width for the calculation, (2) the fraction of rest mass energy that star 
particles return to the gas, (3) the redshift range over which the UV 
background is introduced, (4) whether the abundance of or cooling from H2 
is tracked, (5) and the maximum effective metallicity for metal cooling. 



We also apply a simple self-shielding correction to the 
Lyman- Werner background as described in Shang, Bryan, & 
|Haiman| (pQlO): this uses the local H2 densities multiplied 
by the local Jeans length to determine an estimate of the 
column density of H2 that is shielding each cell. The simple 
estimate in |Draine &; Bertoldi ( 1996 ) is then used to compute 
a self-shielding correction factor. [Wolcott- G reen, Haiman, | 



dynamical times as detailed in Tasker & Bryan ( 2006 ) . The 



Sz Bryan] ( |2011 ) have recently shown that this can result 



in an overestimate of the self-shielding correction for our 
conditions; however, given the highly uncertain level of the 
Lyman- Werner background, this approximation is sufficient 
for this work. 



2.4 Star formation and supernova feedback 

Star particles are created during the course of a simulation 
following the prescription of Cen & Ostriker ( 1992). The im- 



plementation in Enzo is described in more detail in [Tasker 
& Bryan (2006), but we briefly outline it here. Star particles 



are created in a cell if its density rises above 10 5 times the 
universal mean density at that redshift, if the divergence of 
the local velocity is negative, and if the local cooling time is 
less than the dynamical time (or the gas temperature is less 
than 10 4 K). If these requirements are met, then the star for- 
mation rate is computed according to psfr = esFPgas/tdyn, 
where esF — 0.01 is t he star formation efficiency (e.g. 



Krumholz & Tan 



(2007)) and t dyn = (3tt/32Gp) 1/2 is the 



dynamical time. A star particle is then created with mass 
M* = psfrA£Ax 3 , where At is the timestep. To prevent 
very small star particle masses, a stochastic prescription is 
used such that star particles are not created until 100 Mq 
of stars are predicted to form; once that limit is reached, a 
star particle of 100 M© is created (up to a maximum of 80% 
of the baryonic mass in the cell). 

Once they are created, star particles return both en- 
ergy and mass to their surrounding gas, also following the 



prescription of Cen & Ostriker (1992). Star particles return 



some fraction of their rest mass energy to thermal energy in 
the gas. We chose this value to be 3.7 x 10 ~ 6 for most of our 
simulations. This corresponds to 150 M© worth of stars to 
produce one supernova's worth of energy (10 51 ergs). This 
energy is injected as thermal energy into the gas over a few 



energy and mass return fractions depend sensitively on as- 
sumptions regarding the IMF; the values adopted in this pa- 



per are in the ranges determined by other works (e.g. Jung- 



wiert, Combes, & Palous 2001; Kroupa 2007; Oppenheimer 
& Dave|2008| |Conroy fe Gunn|2010| ). Note that we do vary 
the energy feedback return rate in order to investigate the 
impact this has on our results. The particle also returns 25% 
of its mass to the gas and adds 2% of its mass in the form 
of metals. 

Our choice of a purely thermal feedback model is moti- 
vated by simple models of individual supernovae which show 
that the supernova driven sedov blast wave reaches a radius 
of a few parsecs in dense gas before radiative cooling begins 
to slow it down and radiates away energy ( | Che valier|| 1 9 74 \ . 
One might expect that in a galaxy simulation with a resolu- 
tion comparable to this scale, artificial radiative losses seen 



in previous lower resolution studies (e.g. Katz, Hernquist, & 
Weinberg 1992) would be small. In such a simulation this 
resolution would negate the need for artificial subgrid mod- 
els necessary to compensate for this effect (e.g. turning off 
radiative cooling). Our main simulations have a spatial res- 
olution of 10.8 comoving pc; therefore, for z ^ 3, our sim- 
ulations have a spatial resolution less than 3 pc, which is 
the level of resolution needed to test this hypothesis. Given 
the observed old stellar populations found many dSphs and 
their small sizes, these systems are an ideal test case for this 
type of model run at these resolutions. We will discuss the 
efficacy of this model and our assumptions in Section |4.2| 

The 2% metal yield we use is a fiducial value consis- 
tent with the value derived from a Salpeter IMF and mass- 



dependent chemical yields of individual supernovae ( |Madau 
et al.|1996j ) . Uncertainties in individual stellar yields and the 



IMF make this value very difficult to constrain. Moreover, 
more sophisticated chemical model could be employed that 
includes different types of supernovae with different injection 
timescales and chemical yields. We are interested in prob- 
ing the broad effect of a using a purely thermal supernova 
model at high resolution in an AMR code and understand- 
ing the relative effects of supernovae and reionization. For 
these reasons we have chosen to use a simplistic chemical 
model. 
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2.5 Halo tracker and data analysis 



The problem of finding and tracking inheritance between 
dark matter haloes within cosmological simulations is not 



a trivial one. We used the halo finder hop (Eisenstein & 
Hut|1998 ) to identify haloes in our simulation and a simple 
halo tracker to track inheritance between haloes, starting 
with the final target halo at z = and working backwards 
in time to higher redshifts. For each child halo, the tracker 
finds parent haloes in the preceding output by identifying 
haloes that contain dark matter particles that end up in the 
child halo. Parent haloes that contribute more than 10% of 
the child halo's dark matter mass or contain more than 10 7 

in dark matter mass are considered significant parents. 
The lineages of these significant parents are then tracked 
back to previous outputs in the same manner. 

The ability of this simple halo tracker to separate halo 
lineages is limited by the ability of hop to separate interact- 
ing haloes and by the cadence of simulation outputs. There 
were cases of interacting haloes where hop was unable to 
properly separate haloes and the halo tracker was therefore 
unable to properly separate lineages containing these haloes. 
In these cases, we present the potential lineages of interact- 
ing haloes as equivalent and do not attempt to distinguish 
between them. 

We analysed the baryon quantities of identified haloes 
with the simulation data analysis and visualization tool yt 
( Turk et al.|2011 ). We also used yt to construct the images, 
phase diagrams and radial halo profiles which we present. 



3 RESULTS 

We have run a suite of simulations of a single dark mat- 
ter halo, varying the physical prescriptions outlined in Sec- 
tion HI Tabled] summarizes the simulations discussed in this 
paper. Two of our runs, RIO and RIO-earlyUV, include all 
of the physics we have described previously, but differ in 
the time at which the UV backgrounds are introduced. We 
call these two runs our canonical runs, and describe them 
in detail below. We have also run a simulation completely 
neglecting the UV backgrounds, simulation RIO-noUV. This 
simulation, unlike our other simulations which were run to 
z = 0, was run only to z = 3.1. We chose to end RIO-noUV 
at this redshift because of the copious cold, dense gas pro- 
duced. This gas made the simulation run quite slowly and 
by redshift 3.1, the difference in evolution was clear. Sim- 
ulations RIO-lowFB, RIO-noFB and RIO-noFB-LimCool all 
examine the impact of changing the strength of supernova 
feedback. We also tested the importance of H2 for cooling by 
removing it as a coolant in simulation R10-noH2. This set 
of simulations gives us broad levers on the relative impor- 
tance of most of the heating and cooling effects important 
in dwarf formation. In particular, we aim to understand the 
fundamental question of how the dual effects of supernova 
feedback and reionization shape the baryon fractions and 
star formation rates of dwarf galaxies. 



3.1 Canonical runs 

In this section, we describe in detail the results of RIO 
and RIO-earlyUV, which are the highest resolution simu- 
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Figure 1. Top: Evolution of the dark matter mass within r2oo 
(top) and virial temperature (bottom) of progenitor haloes in sim- 
ulation RIO. Coloured lines indicate haloes where star particles 
formed, while gray lines indicate haloes that remain dark. A line 
begins when a halo becomes massive enough to be detected by 
our halo finder and ends when the halo merges into a more mas- 
sive halo. We note that substructure can occasionally separate 
far enough from its parent halo to be detected for a brief time as 
a separate halo, and therefore appears as short lines. In particu- 
lar, two star forming progenitors have a series of close encounters 
during which our halo finder was unable to distinguish between 
them - their evolution is shown in orange. 



lations to include all of the physics outlined in Section [2] 
RIO-earlyUV differs from RIO in that the uniform UV back- 
grounds were turned on in the same way as described in Sec- 
tion |2] but between redshifts 8 and 8.9 instead of between 
redshifts 6 and 7. The purpose of introducing the global 
UV background at different times is to explore the effect of 
patchy reionization. More isolated regions of the universe 
farther from major sources of ionizing photons may be af- 
fected by the ionizing background later than less isolated 
regions. In both RIO and RIO-earlyUV, the universe is ion- 
ized by redshift six. 
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Figure 2. Evolution of baryon properties, within r2oo, of progenitor haloes in simulation RIO (left column) and RIO-earlyUV (right 
column). Top row: The ratio of baryon mass to total mass. Middle row: The ratio of gas mass to total mass. Bottom row: The total 
stellar mass. The cumulative mass distribution of star particles that end the simulation in the final halo is shown in black. The solid line 
colours have the same meaning as in Figure^ The dotted lines bracket the redshift range over which the photoionization background is 
introduced and ramped up to full strength. 



3.1.1 Glo bal properties 

The halo we have chosen is fairly isolated at z — 0; however, 
like all dark matter haloes in cosmological simulations, it 
assembles hierarchically. Figure^ shows the evolution of the 



dark matter halo masses and virial temperatures (as denned 
in Machacek, Bryan, &; Abel 2001) of progenitor haloes in 



RIO and RIO-earlyUV. We present the dark matter evolu- 
tion only for haloes in RIO since the dark matter evolution 
is virtually identical in RIO-earlyUV. At z — 9, there are 30 
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Figure 3. Evolution of the peak density, metallicity and star- format ion rate within progenitor haloes in RIO (left column) and R10- 
earlyUV (right column). Top row: Evolution of the maximum density within each progenitor halo. The density threshold for star formation 
is shown as a dashed line. Middle row: Evolution of the metallicity of the peak density cell shown in the top row. Bottom row: The 
star-formation rate within the last 22 Myr for each progenitor halo is plotted with colour lines, while the star-formation rate inferred 
from the star particles that end the simulation in the final halo is shown in black. The dotted lines bracket the redshift range over which 
the photoionization background is introduced and ramped up to full strength. 
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progenitor haloes more massive than 10 6 M©. Two of these 
haloes are more massive than 10 7 M©, and these two haloes 
gradually build up their mass within two groups over the 
course of the simulation. The two haloes merge at z — 1.8 
in a merger that is about 2:1 in dark matter and nearly 1:1 
in stellar mass in both RIO and RIO-earlyUV. 

The evolution of a variety of baryon quantities in pro- 
genitor haloes is shown in Figures [2] andJ3] We track gross 
properties of progenitor haloes in Figure [2]such as the total 
baryon fraction, gas fraction and stellar mass and quantities 
associated with the densest cell in each halo in Figure [3] such 
as its density and metallicity. 

The evolution of the gas fraction in progenitor haloes 
appears to be dominated by reionization. Figure [2] shows 
sharp declines in the gas and baryon fractions in both RIO 
and RIO-earlyUV at their respective times of reionization. 
These declines are due to photo-evaporative outflows trig- 



gered by reionization (Barkana & Loeb 



1999 



Gnedin h 



Kravtsov 2006). We see that once the gas fraction declines 



during reionization, it remains suppressed for the remainder 
of the simulation. We see no evidence for re-accretion of gas 
once the main halo has been assembled at z — 1.8. 

There are also smaller, but still significant, declines in 
the gas fraction prior to reionization in several luminous pro- 
genitors (Figure |2]). These declines appear to be correlated 
with peaks in the star formation rate as shown in the bot- 
tom panels of Figure [3] Therefore, we conclude that they 
are primarily due to supernova feedback. 

The gas fraction indicates how the bulk of the gas re- 
sponds, but star formation depends on the presence of dense 
gas. The evolution of the density of the densest cell in each 
progenitor halo shares some characteristics with the evolu- 
tion of the gas fraction, but it is not identical. The peak cell 
density declines after reionization, but the decline is delayed 
relative to the declines seen in gas fraction. This is due, in 
part, to self-shielding from the UV background. In RIO par- 
ticularly, there appears to be an extended epoch of elevated 
gas density in one progenitor that does not end until z = 4 
and another progenitor is able to briefly re-cool some gas, 
although this dissipates by z = 3.5. We see that the evolu- 
tion of the peak cell density tracks the star formation rate 
in individual haloes, also shown in Figure [3] which is unsur- 
prising as the presence of dense gas controls star formation 
in our model. 

There are multiple star forming progenitor haloes in 
both simulations; however, RIO-earlyUV has fewer, due to 
the earlier onset of reionization. It appears that some haloes 
which cooled later in RIO do not have time to cool before 
reionization occurs in RIO-earlyUV. The bottom panels of 
Figure [3] show both the star formation rate inferred from 
the stellar population in our target halo at redshift zero, 
and the instantaneous star formation rates in all progenitor 
haloes measured at each data output. The star- formation 
rate inferred from the stellar population at redshift zero is 
built up from star-formation events in multiple progenitor 
haloes with varying durations and intensities. 



3.1.2 A detailed look at the evolutionary stages 

In order to get a better qualitative understanding of how our 
simulated dwarf galaxy evolves, we next look at the spatial 
properties of the halo at various times in RIO. The evolu- 



tion of gross halo properties presented in Figures [T||31 can be 
tracked through the progression of projections presented in 
Figure [4] which are centred on one of the two massive, star- 
forming progenitor haloes of our final halo. We can see the 
spatial variations in gas density, temperature and metallicity 
produced by the processes of cooling, supernova feedback, 
reionization and self-shielding. We can also see the hierar- 
chal build-up of both dark matter mass and stellar mass 
which constitute the final halo. We can break down the evo- 
lution of this halo into a number of states described by these 
projections and by the density-temperature phase plots for 
the same progenitor halo presented in Figure [5] 

The first stage presented is at z = 9 and is shown in 
the top left panel of Figure [5] We see the halo just prior 
to the onset of star formation. Near the virial radius the 
gas shock heats and creates a reservoir of gas at roughly 
the virial temperature of the halo. This gas can then begin 
to cool via H2 cooling (e.g. Abel, Bryan, & Norman 2002). 
Figure [I] shows the evolution of the virial temperature of 
many progenitor haloes; prior to reionization in both the 
RIO and RIO-earlyUV runs, all progenitor haloes have virial 
temperatures less than 10 4 K. Therefore, for primordial gas 
in these haloes, the only important coolant is molecular hy- 
drogen; atomic line cooling is not significant. This stage is 
also shown in projections in the first row of Figure [4] 

Eventually the gas reaches the second stage (top-right 
panel of Figure [5] z = 7.88) where gas temperatures in the 
centre fall to 200 K and gas densities cross 10 -21 g cm -3 , 
roughly the gas density threshold for star formation in our 
model at high redshift. Once gas becomes dense enough to 
form stars in this second stage, we see evidence of supernova- 
heated gas with temperatures in excess of 10 6 K and the 
picture becomes more complex. 

In the third stage (see panel labelled z = 7.04 in Fig- 
ure [5j , gas continues to be heated by supernovae and shows 
evidence of pressure equilibrium at low densities and high 
temperatures, as shown in the phase diagram. Also, with 
the injection of metals, gas can now cool to temperatures 
below 100 K at a broad range of densities. In some cases, 
the gas does appear to be significantly affected by supernova 
feedback in this stage. Figures [2] and [3] show the evolution 
of the gas fraction and the peak gas density of many haloes. 
Prior to reionization there are examples of decrements in 
both quantities, which must be due to supernova feedback. 
The second row in Figure^] also shows this stage. We can see 
that, by this point, metals have been ejected from progen- 
itor haloes and have significantly enriched the intra-group 
medium. 

After z — 7 (in the R10 run), the onset of photoioniz- 
ing and photodissociating backgrounds causes a dramatic 
change in the phase structure of gas within progenitor 
haloes. In this fourth stage (see panels labelled z = 6.36 and 
z = 5.82 in Figure [5), we see nearly all of the gas heated 
to a temperature of about 10 4 K. In some haloes (including 
the one shown in Figures [4] and [5]) , the gas is able to remain 
dense or briefly recool, producing an additional wave of star 
formation. The third row of Figure ^] shows this stage. This 
stage does not appear in haloes in RIO-earlyUV and it ends 
in all progenitor haloes in R10 by z = 3.5. 

After this point, once the dense star- forming gas has 
been dissipated, a small amount of high-temperature gas 
persists within the supernova heated phase in a fifth stage 
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Figure 4. Thin projections centered on a massive star-forming progenitor halo at z = 9.0, 7.04, 5.82, 3.48, and the final halo at z = 
(going from top to bottom). From left to right, the panels show projected gas density, temperature, metallicity, dark matter density, and 
stellar density. The first four redshifts shown have widths of 15 kpc comoving, while the projections at z = have a width of 2r2oo — 
48 kpc. All the projections have a depth of 5 kpc comoving. Projections of gas density, temperature and metallicity are weighted by gas 
density and projections of dark matter density and stellar mass density are unweighted. 



(see panel labelled z = 4.97 in Figure [5]). However, this gas 
phase eventually dissipates and the small amount of remain- 
ing gas sits at a temperature slightly above 10 4 K. This is 
both the temperature of the IGM and roughly the virial 
temperature of the final halo. 



At low redshift (see the last two panels of Figure [5]), a 
series of dry mergers continue until the final halo is assem- 
bled. The final haloes that result in RIO and RIO-earlyUV 
are gas poor spheroids of dark matter and star particles. RIO 
and RIO-earlyUV differ in stellar mass by an order of magni- 
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Figure 5. Density-temperature distributions of gas within the r2oo radius for one of the massive halo progenitors at the redshifts shown. 
The colour of cells shows the cumulative sum of mass within the corresponding density-temperature bin. 



tude: the halo in R10, which had a late onset of reionization, 
has a final stellar mass of 1.43 x 10 6 M©, while the halo in 
RIO-earlyUV, which had an earlier onset of reionization, has 
a final stellar mass of 1.16 x 10 5 M . Table [| summarizes 
some of the final halo properties in these simulations. 



3.2 Alternate physics runs 

We have also conducted a number of simulations in which 
we have altered some aspect of our physical prescriptions in 
order to better understand their effect on the evolution of the 
final halo. We have done simulations that have changed some 
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Table 2. Summary of Final Halo Properties. 





RIO 


RIO-earlyUV 


fl/T , , / A/f^ 
lV-ltot/-LVl@ 


1.00 X 1U 


1.00 X 1U 


M*/M 


1.43 x 10 6 


1.16 x 10 5 


MgaslMtot 


3.24 x 1CT 3 


3.25 x 10- 3 


r 2 oo (kpc) 


23.7 


23.9 


^1/2 (PC) 


704 


213 


M 1/2 /M 


3.05 x 10 7 


3.86 x 10 6 


M 3O o/M 


7.53 x 10 6 


7.41 x 10 6 


cr 1/2 (km/s) 


7.83 


8.30 


Log((Z/Z >) (median) 


-0.29 


-1.22 


Log((Z/Z » (mean) 


-0.076 


-0.92 




0.84 


0.14 



Note: The quantities presented in each row are (1) the total 
mass within r2oo, (2) the total stellar mass within r2oo, (3) the 
total gas fraction within r2oo, (3) r200, the radius within which 
the mean halo density is 200 times the critical density of the 
universe, (4) the radius enclosing half the stellar mass, (5) the 
total mass within 7*1/2? (6) the total mass within 300 pc, (7) the 

velocity dispersion of star particles within 7*i/ 2 , (8) the 
mass-weighted median of the star particle metallicities, (9) the 
mass- weighted mean of the star particle metallicities, (10) the 
mass-weighted standard deviation of the star particle 
metallicities. 

aspect of the cooling prescription, the supernova feedback 
prescription and the UV background prescription. Table [I] 
summarizes these simulations. 



3.2.1 Simulation without iJ 2 cooling 

We conducted a simulation without the creation of or cool- 
ing from H2. Simulation R10-noH2 has line cooling from 
six atomic species of hydrogen and helium as well as metal 
line cooling and heating from the Cloudy models of |Smith^] 



Sigurdsson, fe Abelj ( [2008] ) . The UV backgrounds, the self- 
shielding by neutral gas, and the star formation and feed- 
back models are the same as RIO. This simulation produced 
an entirely dark halo; no star formation occurred. We found 
that without molecular hydrogen cooling, gas temperatures 
in the most massive progenitors did not fall below 500 K 
and densities remained below 6 x 10 -24 g cm -3 . The peak 
in gas density occurred at z = 5.8, shortly after reionization 
reached full strength. Our model requires gas to be denser 
than 8.09 x 10 - 23 g cm -3 to form star particles at this red- 
shift. With no initial star formation, there is no source for 
metals, which are the most important gas cooling agents for 
star formation. Once the photoionizing background is turned 
on, the halo gas is further heated and the gas becomes even 
more difficult to cool. The result is an entirely dark halo at 
z = 0. 



gross differences in the properties of the final stellar compo- 
nent between these two simulations. The total stellar mass 
in the final halo differs by an order of magnitude between 
the two runs and the star-formation history is more bursty 
and truncated in RIO-earlyUV. The median stellar metallic- 
ity is much lower in RIO-earlyUV (0.06 Z versus 0.51 Z ), a 
consequence of the lack of late and extended star-formation 
bursts. 

To further explore the effect of the UV heating on the 
halo's evolution, we conducted a simulation without the pho- 
toionizing and photodissociating UV fields, RIO-noUV. This 
simulation was only run to a redshift of 3.1, but interesting 
comparisons can still be made. Figure [6] shows the evolu- 
tion of a variety of gas properties for the two massive pro- 
genitors of our final halo in RIO-noUV (we show only two 
haloes for clarity). There are some striking similarities and 
differences in the evolution of these haloes between R10 and 
RIO-noUV. We see dips in the gas fractions of haloes in R10- 
noUV, some quite dramatic, but the gas fraction is able to 
rebound in both massive progenitors. By z = 3.1, when we 
stop RIO-noUV, there is an order of magnitude difference in 
baryon and gas fractions for both haloes between the two 
simulations. The baryon fraction in one halo has rebounded 
to its level at z = 12. 

The dense gas exhibits similar behavior. The peak gas 
density does drop considerably between redshifts of 6 and 
4, but the gas is able to re-cool and condense to a level 
where it is able to form stars. This produces an uptick in 
the star formation rate and several haloes are forming stars 
at z = 3.1 when the simulation was ended. 

The dips seen in the gas fraction and peak gas density 
in RIO-noUV are most likely due to the effect of supernova 
feedback. In the absence of UV heating, supernovae act as 
the main regulator of the star formation rate in this run. The 
degree of supernova heating scales with the star formation 
rate and once the rate drops, gas is able to both re- accrete 
onto the halo and re-cool in its centre. This produces an 
increase in the star formation rate. We did not run this 
simulation beyond z = 3.1, so it is impossible to say whether 
this cyclical process will continue or whether the constant 
bursts of star formation will produce an overall net decline 
in the gas fraction or peak gas density. 

These simulations indicate that the UV background 
plays a crucial role in regulating the gas fractions of the 
progenitor haloes. In our model, it produces a sharp drop 
in the gas fraction in all haloes. This drop leaves an imprint 
on the star formation history inferred from the final stellar 
population of our target halo (the solid black lines in the 
bottom two panels of Figure [3| . This imprint is particularly 
evident because our UV backgrounds are spatially uniform, 
so all progenitors are affected with the same level of flux, 
and all the progenitor haloes are in the mass range sensitive 
to the UV backgrounds. 



3.2.2 Simulations with alternate UV backgrounds 

Our canonical simulations demonstrate that the UV back- 
ground has an important effect on the evolution of the halo 
gas. We observe that the gas fractions of progenitor haloes 
drop dramatically at the redshift of full-strength reioniza- 
tion in both our early reionization simulation (RIO-earlyUV) 
and our late reionization simulation (R10). We also observe 



3.2.3 Simulations with alternate supernova feedback levels 

We also conducted three simulations exploring the effect of 
changing the strength of supernova feedback in our model. 
In most of our simulations, including R10, star particles re- 
turn 3.7 x 10 ~ 6 of their rest mass energy to thermal gas en- 
ergy. In simulation RIO-lowFB, we have set this fraction to 
be 10~ 6 . In simulations RIO-noFB and RIO-noFB-LimCool, 
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Figure 6. Evolution of properties for two massive progenitors in simulations RIO (black) and RIO-noUV (green). RIO-noUV was not 
run past a redshift of 3.1. The properties shown here are the halo baryon fraction within the r2oo radius (top left); the halo gas fraction 
within the r2oo radius (top right); the density of the densest cell within each halo (bottom left); and the stellar mass within the r2oo 
radius for each halo (bottom right). The density threshold for star formation is shown as a dashed line in the bottom left panel showing 
the peak gas density. 



star particles return zero thermal energy to the gas through 
supernovae, however, star particles do return metals to the 
gas in the same way as in RIO and RIO-lowFB. 

Properties for the two most massive progenitor haloes 
in simulation RIO-lowFB are shown in Figure [7] The gas 
fractions for these two haloes behave in nearly the same 
fashion as they do in simulation RIO. There appear to be 
smaller declines in the gas fraction prior to z = 6, which 
must be due to the reduced impact of supernova heating. 
Curiously, while the gas fraction declines in the same way 
and to the same levels as in RIO, the peak gas density within 
each halo behaves differently. There are sharp declines in 
the peak gas density between redshifts 6 and 4, but the 
peak gas density is then able to rise significantly in one of 
the progenitors, up to the level required for star formation. 
There appears to be a steady level of star formation in R10- 
lowFB down to low redshift, unlike in RIO. 

Star particles in our model impact the thermal state 
of the gas in two ways: first, by heating it through super- 



nova feedback, and second, by injecting metals that act as 
coolants. In simulations RIO-noFB and RIO-noFB-LimCool 
we have completely eliminated thermal feedback while main- 
taining the metal feedback. Unsurprisingly, we see metals 
collecting in the centres of haloes to a much greater degree 
than in RIO, making it unclear if changes in the no-feedback 
runs are due to a lack of energetic input, or an increase in 
the metal cooling rate. 

To differentiate between these causes, in run RlO-noFB- 
LimCool we capped the effective metallicity used for cooling 
at 0.1 Z . Figure [8] shows that the differences between R10- 
noFB and RIO-noFB-LimCool are minimal; therefore the 
differences seen between these two simulations and RIO are 
due to a lack of heating by supernovae, not to enhanced 
cooling from the collected metals. 

The gas fractions of haloes in RIO-noFB and RlO-noFB- 
LimCool decline in almost the same way as in RIO after 
the introduction of the UV fields. They are slightly elevated 
down to z ~ 2, as compared to the gas fractions in RIO, 



© 2012 RAS, MNRAS 000, [l]-?? 



14 C. M. Simpson et al 




6 8 
Redshift 



10" 



c 
o 



u 



10' 



10 _; 




6 8 
Redshift 




4 6 8 10 12 
Redshift 




6 8 
Redshift 



Figure 7. Evolution of properties for two massive progenitors in simulations RIO (black) and RIO-lowFB (blue). The properties depicted 
are the same as those shown in Figure [6] 



but this may be due to their higher gas fractions at z = 6 
when the ionizing background reaches its full strength. As 
in RIO-lowFB, the overall decline in the gas fraction does 
not correlate with an overall decline in the peak gas density. 
The peak gas density in both massive progenitors dips at 
z = 6, but it does not fall below the threshold density for 
star formation, and we see that a low level of star formation 
persists in both progenitors down to z = 0. 



In all four simulations that explore different levels of 
supernova feedback (RIO, RIO-lowFB, RIO-noFB and R10- 
noFB-LimCool), the final gas fraction is extremely low and 
virtually identical. Comparing these results with those from 
RIO-noUV, we can conclude that the factor that mainly de- 
termines the gas fraction of the halo is the UV background. 
Supernova feedback does play an important role in regulat- 
ing the dense gas component, from which stars are formed. 
UV heating can produce an initial, sharp suppression of the 
peak gas density, however it is supernova heating that ap- 
pears to determine whether it remains suppressed. 



4 DISCUSSION 

We have explored a variety of physical effects in a series of 
high-resolution, hydrodynamical, cosmological simulations. 
Our goals were to better understand the physics that regu- 
late star formation in a dwarf halo comparable in mass to the 
dwarf satellite galaxies seen in the Local Group, as well as 
to test current feedback models in such systems. We found 
that the combination of supernova feedback and reioniza- 
tion causes a large suppression of the star formation rates 
in these haloes. Reionization appears to regulate the overall 
gas budget within the halo. Supernova feedback seems to 
more directly affect the dense gas that is available for star 
formation. Of course, these heating sources are not strictly 
exclusive in their effects; supernova feedback can have a 
quasi-periodic effect on the gas fraction and reionization can 
cause a dip in the peak gas density. However, the quantities 
affected are typically able to recover in the absence of their 
primary regulator. 

The reason for this dichotomy between reionization and 
supernova feedback is primarily their different physical loca- 
tions within the halo. The UV backgrounds we implement 
are uniform and are attenuated by our self-shielding pre- 
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Figure 8. Evolution of properties for two massive progenitors in simulations R10 (black), RlO-noFB (yellow) and RlO-noFB- 
LimCool(red). The properties depicted are the same as those shown in Figure [6] 



script ions in dense gas at halo centres. At the onset of reion- 
ization, the ionizing background is able to produce a pressure 
imbalance nearly everywhere in the halo except the very cen- 
tre. A large fraction of gas is therefore driven out of the halo 
in photoevaporative winds, while the central dense gas can 
remain intact for some time (Alvarez, Bromm, & Shapiro 



2006]|Gnedin fe Kravtsov||2006| jWhalen et al.|2008 ) 



Supernova feedback only injects thermal energy into 
cells containing young star particles in our model. Since 
the star particles tend to form and remain centrally located 
in haloes throughout the simulation, the energy they inject 
tends to affect central dense gas. We do see evidence that 
supernovae in our model affect regions outside the central 
dense gas in the form of supernovae driven outflows which 
enrich the IGM with metals, as seen in Figure [4] and dis- 
cussed more below. It is unclear, however, to what extent 
our supernova feedback model accurately captures the heat- 
ing of the ISM. While the simulations that we have run 
with little or no thermal feedback seem to indicate that our 
feedback prescription is key in expelling dense gas and trun- 
cating the star formation history, we will explore this issue 
in more detail later in this section. 



In the following sections, we will first compare our re- 
sults to observed properties of low mass dwarf systems, then 
we will examine the effect of supernova heating in more de- 
tail, followed by discussions of spatial resolution and metal 
dispersal, and ending with a discussion dark matter proper- 
ties and neglected physics. 



4.1 Observable properties 

In this section we compare our models to observations of 
several types of low- luminosity dwarf galaxies. This com- 
parison is important because understanding the ways in 
which our simulations fail to produce realistic galaxies can 
provide important insights into the specific ways in which 
our model prescriptions fail. We can consider two classes of 
dwarf galaxies for comparison: mass analogues and environ- 
mental analogues. Our simulated halos have stellar masses 
and velocity dispersions similar to those of dSphs found near 
the Milky Way. Unfortunately, their proximity to a large 
galaxy does not make them suitable environmental ana- 
logues; our simulated halo is 2.9 Mpc from a Milky Way 
sized halo. Dwarf galaxies have been observed in underdense 
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Figure 9. Comparison between the observational properties of 
Local Group dSphs and the simulated properties of the final 
haloes produced in R10 and RIO-earlyUV. All quantities are plot- 
ted against total V band luminosity. V band luminosities are esti- 
mated for our simulated haloes as described in section [4. Top 
panel : Velocity dispersio ns for 28 Local Group dSphs collected 
by |Walker et al.[ ( 2009b but including the velocity dispersion 
for the dynamically cold component of Bootes I measured by |Ko-| 
|posov et al.|(20li} (open circles); and the three-dimensional veloc- 
ity dispersions of star particles within the half-stellar mass radii 
of R10 and RIO-earlyUV (solid squares). Middle panel: The half- 
light radii for the same sample of dSphs (open circles) and the 
half stellar mass radii for the final haloes in R10 and RIO-earlyUV 
(solid squares). Bottom panel: The mean stellar metallicities for 
15 Milky W ay dSphs measured by |Kirby et al. 1(2008) and |Kirby| 
|et al.|(2011) (open circles), the mean (solid squares) and median 
(solid triangles) star particle metallicities for the final haloes in 
R10 and RIO-earlyUV. 



environments, however, most of these are dlrrs that are more 
luminous than both the typical Milky Way dSphs and our 



simulated haloes (Weisz et al. 2011a). A handful of dSphs 



are known in low-density environments (e.g. Makarov et al. 



2012] ), however while some of these systems may be relics 

others likely under- 



of reionization (Monelli et al. 



2010) 



went pericentric passages with a large host galaxy ( Teyssier, 



Johnston & Kuhlen 2012). Many of the lowest luminosity 



dwarfs discovered around the Milky Way have been found in 
the SDSS, for which the distance limit for detection has been 
demonstrated to be less than the Milky Way virial radius 
(Koposov et al. 2008). It is therefore possible that similar 



low- luminosity dwarfs await detection at greater distances 
in low-density environments. 

Our canonical simulations R10 and RIO-earlyUV pro- 
duce a final halo similar in some respects dSphs of the Local 
Group, although in other respects our models differ. Table|2] 
lists some of the properties of these haloes. The final velocity 
dispersion of star particles within the radius enclosing half 
the stellar mass, r 1 / 2 , is roughly 8 km s _1 in both simula- 
tions. This makes them dynamical analogues of Milky Way 
dSphs, which have velocity dispersions ranging from 3 to 16 
km s _1 as shown in Figure [9] Both simulated haloes have 
masses within 300 pc of approximately 7.5 x 10 6 M . This 
value is close to 10 7 M©, which is the typical dynamical mass 
derived for Milky Way dSphs from velocity measurements of 
stars within 300 pc ( |Strigari et aL|2008 ). 



4-1.1 Total stellar mass 

The final stellar masses produced in our simulated halo 
are 1.43 x 10 6 M in R10 and 1.16 x 10 5 M in R10- 
earlyUV. Assuming a stellar mass-to-light ratio in the V 
band of 2 M /L for our (old) stellar population (Krui- 



jssen fe M ieske 2009), the haloes in R10 and RIO-earlyUV 
would have V band luminosities of approximately 7.15 x 10 5 
L and 5.8 x 10 4 L respectively. These estimates place 
the halo in RIO-earlyUV at the upper luminosity end of the 
ultrafaint dSphs and R10 among the lower luminosity clas- 
sical dSphs (Walker et al]|2009b| [Simon fe Geha||2007| ), as 
shown in Figure [9] If we compare our estimated luminosi- 
ties to the haloes' total gravitational masses, these haloes 
have mass-to- light ratios within 7*1/2 (the radius enclosing 
half the stellar mass) of 85 M /L for R10 and 133 M /L 
for RIO-earlyUV. The spatial extent of the final stellar com- 
ponent is also comparable to real systems; Figure [9] shows 
that the half stellar mass radii for the final haloes in R10 
and RIO-earlyUV are similar to the half-light radii of dSphs 
with similar luminosities. 

In other highly resolved calculations of low-mass dwarf 
haloes, the total stellar mass is generally similar to what 



Hatzidimitriou 1 (1995) 
Walker et ah] (2009c); | Mart in, de Jong, fe Rix| (2008); | Walker 



These data were ta ken by |lrwin 



et al.N2007^|Mateo, Olszewski, fe Walker| (2008); |Koch et al. 



(2007 



(2007 



Martin et al.| (2007); |Koch et a l. ( 2009 ); |Simon fe Geha 
Aden et al. (2009); B elokurov e t al. (200 8); | Walker et al 



( 2009a) ; |Irwin et al.j (2007); |Geha et al . (2009 ); Belok urov et al. 



( 2009) ; |McConnachie !xwin ( 2006) ; |C6te et al.| (1999) ; [Chap- 
man et al.N2005);|Ibata et al.M2007);|Letarte et al.|(|2009);| Lewis 
et al. 1(2007) ;|Ibata et al.|(2007);|Majewski et al.|(2003) ; |Saviane. 
Held, fe Piotto|(1996) ; |Fraternali et al.|(2009) 
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we have found. [Sawala et al. (2010) found, for their highest 
resolution model of a similar halo mass, a simulated stellar 
mass of 7.35 x 10 6 M©. Recent work by Governato et al 



(2012) also probed halo masses close to 10 M©, and while 
the stellar content is not the focus of their study, their mod- 
els appear to produce only ~ 5 x 10 4 M© of star particles for 
haloes close to this mass. One of the many differences be- 



et al. 


(2010) 


et al. 


(2012} 



reionize the universe at z ~ 6 and Governato 
do so much earlier at z ~ 9. Considering this, 
the results of | Sawala et al.| ( |2010| are most analogous to our 
simulation R10, which has a final stellar mass of 1.43 x 10 6 
M© , and the results of Governato et al. ( 2012 ) are most anal- 



ogous to simulation RIO-earlyUV which has a final stellar 
mass of 1.16 x 10 5 M©. In this light, these results appear to 
be roughly consistent with each other, although differences 
in final stellar mass between analogous haloes can approach 
an order of magnitude. Also, while the simulated haloes are 
similar, they are not identical realizations, and other dif- 
ferences in the physical models being tested (i.e. feedback 
prescriptions) should not be underestimated in their effect. 



4.1.2 Metallicity 

Milky Way dwarfs are observed to lie on a fairly tight 
luminosity-metallicity relation (Kir by et al.|201l] ), shown in 
Figure [9} which may extend beyond the Milky Way ( Kirby, 
Cohen, &; Bellazzini 2012). Figure [9] shows the metallicity- 
luminosity relation found for Milky Way dSphs. Star par- 
ticles in R10 have a mean metallicity of 0.84 Z© and a 
median metallicity of 0.51 Z©; in RIO-earlyUV they have 
a mean metallicity of 0.12 Z© and a median metallicity of 
0.06 Z©. If we make the estimate t hat [Fe/H] ~ log(Z/Z©) 
Woo, Courteau, & Dekel 2008}, then these values trans- 
to mean and median [Fe/H] values of -0.076 and - 
0.29 dex for R10 and mean and median [Fe/H] values of 
-0.92 and -1.2 dex for RIO-earlyUV. Our assumption that 
[Fe/H] no log(Z/Z©) is a somewhat flawed one given the 



(e.g. 
late 



short duration of star formation in our progenitor haloes. 
We would expect this short duration to result in elevated 
values of [a/Fe] for their stellar populations ( Tinsley|T9 79 ) . 
This has the potential to decrease our estimated [Fe/H] by a 
few fractions of a dex, but we are interested in understand- 
ing differences in [Fe/H] of a dex or more. 

Based on our estimated luminosities calculated in Sec- 
tion [XTTT] and according to the dSph luminosity-metallicity 
relation of Kirby et al. ( 2011} , star particles in the final halo 
in R10 should have an average [Fe/H] of -1.76 dex (1.7x 10~ 2 
Z©) and -2.09 dex (8.1 x 10" 3 Z©) in RIO-earlyUV. Fig- 
ure [9] shows this discrepancy, which is over 1.5 dex for R10 
and over 1 dex for RIO-earlyUV. This 1-2 dex difference be- 
tween observed systems and our simulations has important 
implications for the utility of the feedback model we use. It 
appears that while our adopted feedback model effectively 
suppresses the formation of dense gas and stars, it does not 
drive sufficient metal ejection from the ISM. 



4-1.3 Gas content and star formation history 

Our simulated halo in R10 and RIO-earlyUV, has a final gas 
fraction that is very suppressed; the total mass in gas is less 
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Figure 10. Virial temperature versus redshift for progenitor 
haloes that are constituted of at least 50% final target halo dark 
matter particles. Haloes plotted in black contain no star particles; 
haloes plotted in cyan have formed star particles within the last 
20 Myrs; and haloes plotted in red contain star particles, but have 
not formed star particles within the last 20 Myrs. 



than 1% of the total halo mass, well below the cosmic baryon 
fraction (see Table |2j). As we have discussed in Section 3.1.2 



this remaining halo gas is very diffuse and highly ionized 
with a temperature of a few factors of 10 4 K (see the z — 
panel of Figure [5]) and would not be directly observable in 
a real system. 

Many low- luminosity dSphs in the Milky Way are very 
gas poor with little or no H I, but some of the more distant 
Milky Way dwarfs can be quite H I rich, for example Leo 
T ( |Ryan- Weber et aT||2 008). Our simulated halo is isolated 
and has had no interactions with the hot halo of a massive 
galaxy during its formation. We might, therefore, expect it 
to be more analogous in its gas content to the more distant 
Local Group dwarfs. The observed strong correlation be- 
tween environment and gas content for Local Group dwarfs 
(Grcevich &; Putman|[2009 ) suggests that their gas content 



may be due to environmental effects, which this study does 
not explore since we only simulate one environment. 

Reionization appears to be the mechanism that removes 
most of the gas mass in our models and the only simulation 
with evidence for late time gas accretion is RIO-noUV, which 
had no UV background. From this simulation, we conclude 
that the lack of reaccretion of gas onto our simulated halo 
is due to its low virial temperature relative to the warm 
temperature of the IGM after reionization, a consequence of 
the halo's low mass (T V i r oc M 2//3 ). It does not appear that 
supernovae are a factor in preventing the reaccretion of gas, 
as the suppression of gas fractions is remarkably similar in 
our simulations with different levels of supernova feedback. 

Figure [To] shows that at late times the main progenitor 
haloes and massive subhaloes in R10 all have virial tem- 
peratures well below the temperature of the IGM (a few 
factors of 10 4 K). This plot also shows that, prior to reion- 
ization, with a few exceptions, progenitor haloes generally 
need to have virial temperatures above 2000 K in order to 
form stars, due to the inefficiency of H 2 cooling, see Tegmark 
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et al. (1997). After reionization, most progenitor haloes stop 



forming stars. As we discussed in Section 3, this suppres- 
sion of star formation was due to a combination of the large 
scale photoevaporation of gas from progenitor haloes and the 
effect of supernova heating on the remaining self-shielded 
clumps of dense gas. The net effect is consistent with the 
idea that reionization introduces a filtering mass or filtering 
virial temperature, below which star formation is truncated 
(e.g. |Gnedin|2000| ). 

Observations of dwarfs in the present epoch cannot tell 
us the complete evolution of gas in these systems, but their 
star formation histories derived from color-magnitude dia- 
grams do give us the history of dense star-forming gas. We 
find that our simulated star formation histories in RIO and 
RIO-earlyUV are sharply truncated after reionization. RIO 
has the most extended star formation history; however, all 
star formation in this model ends by a redshift of 3.5. Milky 
Way dSphs that are closest in luminosity to RIO and R10- 
earlyUV have a range of velocity dispersions and star forma- 
tion histories. In general, the most truncated star formation 
histories tend to be found in s ystems with the low est veloc- 
ity dispersions (~ 3 km s _1 ) ( |Brown et al.||2012 ), although 
recent analysis suggests at least one system with a lumi- 
nosity and velocity dispersion comparable to our simulated 
haloes, Canes Venatici I, has only ancient stellar popula- 



tions (Okamoto et al. 2012). In general, however, most of 



the Milky Way dwarfs contain some intermediate-age stellar 
population (Tolstoy, Hill, fe Tosi|2009t which we do not see 
in our simulations. 



4-1-4 Implications of observational comparisons 

Since the size of our simulated halo appears to be what 
governs its gas content and star formation history, the lack of 
a known mass or luminosity analogue for our isolated system 
due to the limitations of current observational surveys makes 
evaluating our model's ability to produce systems with the 
proper gas fraction problematic. 

Given these difficulties, what can we conclude about 
the discrepancies between our simulated haloes and observed 
Local Group dwarf systems, especially in regards to the lack 
of an intermediate age stellar population? The star-forming 
state of our simulation appears to be heavily dependent on 
the metagalactic UV background we implement and this is 
relatively well understood for z < 6. The lack of dense, star- 
forming gas in our simulations at late times implies one of 
four things: one, either the implementations of the meta- 
galactic UV background or the local shielding are inade- 
quate and a larger reservoir of gas should survive reioniza- 
tion in these systems; two, our supernova feedback prescrip- 
tion is too effective in destroying dense gas that survives 
reionization; three, the handful of gas rich, but low luminos- 
ity, dwarfs (such as Leo T) observed at large galactic radii 
have virial temperatures above that of the IGM (and above 
the halo we simulate) which allows them to accrete gas at 
late times; or four, our halo merger tree is atypical for this 
halo mass or is atypical for dwarfs that end up in dense 
environments. 

The first scenario is possible and more sophisticated 
methods may be used to simulate ionization and shield- 
ing, as we will discuss in Section |4.6| However, self-shielding 
would have to be many times more effective than in our 



current model to significantly impact the gas fractions. The 
second scenario also seems unlikely, since the discrepancy 
in stellar metallicity between our simulations and observed 
systems indicates that our feedback prescription is not ef- 
fective enough. If the third scenario is true, then we would 
predict undiscovered dwarfs with the luminosities, star for- 
mation histories and gas content similar to Milky Way dSphs 
in the field. The fourth possibility is that our halo merger 
tree is different from dwarfs found in denser environments. 
Although we do not study this directly, this possibility is 
plausible. It is certainly the case that the typical mass of 
progenitor haloes at the time of reionization greatly impacts 
the final properties of the system and this typical mass may 
be different for different merger histories. It is interesting 
to note, however, that analysis of more isolated, but more 
luminous, dwarfs in the local volume indicates general con- 
sistency in star formation histories between these systems 
and Milky Way dwarfs ( |Weisz et aL|2011b ) 



4.2 Supernova feedback model 

Our supernova model involves a simple injection of thermal 
energy into the ISM over a few dynamical times. For our 
canonical models we chose the level of thermal energy to in- 
ject into the gas to be 10 51 ergs per 150 M© of star particles. 
As discussed in Section [5.2.3| we find that lowering the feed- 
back energy, as in simulations RIO-lowFB, produces a more 
extended star formation history with star formation contin- 
uing past z = 1. The suppressed star formation histories in 
RIO and RIO-noUV are generally consistent with other cos- 



mological simulations probing halo masses of 10 9 M©. Gov 



ernato et al. (2012) found star formation to be strongly sup- 



pressed in haloes with viral temperatures below that of the 
IGM and did not see the formation of central dark matter 
cores in haloes of this mass, which were produced in higher 
mass haloes by repeated starbursts extending down to late 



times. iSawala et al. (2010) found star formation lasted ap- 



proximately 2.62 Gyrs in their highest resolution simulation 
closest in mass to our simulated halo, a somewhat longer 
star formation epoch than found in R10. As we have dis- 
cussed, observations of ultrafaint dSphs in the Milky Way 
have found systems with truncated star formation histories 
(Bro wn et al.||2012 Okamoto et al.|[2012 ). However, many 
isolated, but most likely more massive, dwarfs have had star 
formation events since z = 1 ( |Weisz et al.||2011a| ). 

While a lower level of energy conversion may be indi- 
cated by the lack of late time star formation in our sim- 
ulated halo, the metallicity of the star particles produced 
in R10 and RIO-earlyUV exceeds that observed, arguing in 
the opposite direction. To better understand how supernova 
heating occurs in our model, we examine heating and cool- 
ing timescales for each cell. We estimate the heating time 
of a cell to be the ratio of the cell's thermal energy to the 
average supernova heating rate, assuming that exactly one 
star particle is heating that cell, spread out over the typical 
timescale we adopt in the code: 

3 mceiikT ce ii tf ^ 



theat — 



2 HceiirriH e S NM*c 2 ' 



where esN = 3.7 x 10 -6 is the fraction of its rest mass energy 
a star particle returns to the gas; tf is the time over which a 
star particle returns its supernova energy to the ISM in our 
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log Heating Time [s] 

Figure 11. Cooling time versus estimated heating time at z = 
7.04 for cells within r^oo for the massive progenitor whose phase 
evolution is shown in Figure [5] The line where the cooling time 
equals the heating time is shown in black. The colour of cells 
shows the cumulative sum of mass within the corresponding cool- 
ing time-heating time bin. 

model, which we take to be 10 Myr; and M* = 100 M©, a 
typical star particle mass in our simulations. We computed 
theat for every cell within r 2 oo and compare it to the cooling 
time of that cell in Figure [TT] computed at z = 7.04. To 
be clear, the cells for which we have estimated this heating 
time are not necessarily being heated by a star particle at 
the time shown; we are simply estimating a timescale for the 
response to such heating. 

From Figure |11| we see that for nearly all of the cells 
within the halo, the cooling time far exceeds the potential 
heating time. However, we see that gas with the shortest 
cooling times lies along the line where heating time equals 
cooling time. This gas is most likely the star forming gas 
since star particles tend to form in cells with lower cool- 
ing times (more precisely in cells where the cooling time is 
shorter than the dynamical time). Star particles are more 
likely to heat this gas since young star particles do not wan- 
der too far from their formation sites. 

We interpret this result to mean that some energy pro- 
duced by our feedback model is cooled in an unrealistic fash- 
ion. The key issue is most likely the length of time over which 
supernova energy is injected into the ISM, i.e. the length of 
tf. Our feedback algorithm was motivated by the require- 
ments of simulating larger mass halos and was designed to 
model the release of supernova energy from an entire stellar 
population that forms over an extended period of time (more 
than 10 7 years). A typical star particle in our simulations 
has a mass of 100 M©, which is approaching the mass of an 
individual massive star. Extending the release of supernova 
energy over 10 7 years is therefore not suitable for this smaller 



mass star particle which is producing about one supernova's 
worth of energy. Reducing tf or even eliminating it and in- 
stantaneously injecting energy into the cell would be more 
appropriate and has been explored in idealized simulations 
of the ISM ( |Joung & Mac LowpOOGj ). We found, however, 
that even with our high spatial resolution, injecting a super- 
nova's worth of energy into one resolution element instan- 
taneously proved to be computationally challenging and we 
are working to develop an improved algorithm. 

In our current model implementation, energy that is 
able to escape the central regions with shorter cooling times 
will not be quickly lost since gas everywhere outside the cen- 
tral region has comparatively long cooling times. Indeed in 
Figure [4] we do see evidence for supernova heated bubbles 
and metal pollution of the ISM and IGM. Previously in Sec- 
tion [3] we had concluded that reionization determines the 
overall gas fraction of haloes and that dense gas is more di- 
rectly affected by feedback. This interpretation still largely 
holds, however, the effect on dense gas may be more pro- 
nounced than what we have observed in our simulations. 
This would have the effect of enhancing the dominance of 
feedback over the UV background as the regulator of dense 
gas. 

4.3 Resolution effects 

To explore the impact of spatial resolution, we examine sim- 
ulation R43, which is a run with the same dark matter parti- 
cle resolution and physics as R10, but with four times coarser 
spatial resolution. R43 has a minimum cell length of 43 co- 
moving parsecs. Figure ^] shows the evolution of relevant 
quantities for the two main star forming progenitors in R43 
and Table [3] lists many of the final halo's properties. The gas 
fractions of haloes in R43 evolve in a similar way to those in 
R10, with a large decline around z = 6. The declines in gas 
fraction prior to reionization, however, are almost nonexis- 
tent. These declines in R10 were due to supernova feedback 
and we see that coarser spatial resolution has decreased its 
effect. 

The dense gas, which we found was more directly regu- 
lated by supernova feedback, also has very different behav- 
ior. The dense gas in the two main progenitors does not rise 
to quite as high a peak density, but is able to survive longer 
because of the diminished impact of supernovae in the low 
resolution run. Cells in R43 encompass 64 times the volume 
of cells in R10; therefore when a star particle injects energy 
into a cell in R43, the specific heating rate is approximately 
l/64 th of what it is in R10. Dense gas is therefore able to 
survive longer. The maximum density cell within each halo 
peaks at a somewhat higher density in R10 because we are 
able to better resolve the collapse of cold clouds. 

The star formation events in R43 are of longer dura- 
tion due to the enhanced survivability of dense gas. Star 
formation continues in the main massive progenitor down 
to z = 2. This produces a final halo with a greater stellar 
mass (6.82 x 10 6 M©). This is somewhat counter to the ex- 
pectation that increased resolution will produce more star 
formation (since we resolve denser gas), however this result 
makes sense in light of the resolution-dependent effectiveness 
of our supernova feedback model. 

Two other high-resolution, cosmological simulation 
studies track the evolution of similar mass haloes to z = 0: 
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Figure 12. Evolution of properties for two massive progenitors in the canonical simulations R10 (black) and the lower minimum 
resolution simulation R43 (magenta). The properties depicted are the same as those shown in Figure [6] 



Table 3. Summary of Final Halo Properties for Coarse Resolu- 
tion Simulation R43. 



R43 



Mtot/MQ 


1.54 x 10 9 


M*/M© 


6.82 x 10 6 


Mgas/Mtot 


3.28 x 1CT 3 


T200 (kpc) 


23.7 


ri/2 (pc) 


326 


M 1/2 /M 


1.56 x 10 7 


M 3O o/M 


1.42 x 10 7 


cr 1/2 (km/s) 


8.56 


(Z/Zq) (median) 


0.80 


{Z/Zq) (mean) 


1.0 


vz/Zq 


0.83 



Note: The quantities presented in each row are the same as in 
Table 12 



Sawala et al. (2010) and Governato et al. (2012). Both stud- 



ies use a SPH code, so it is not possible to directly compare 
our resolutions in a consistent way. However, both simula- 
tions are closer in spatial resolution to our simulation R43; 



Sawala et al. (2010) fix their softening length in such a way 
that their resolution is typically less than 100 pc and |Gover^| 
nato et al. (2012) have a softening length of 64 pc for their 
highest resolution run. In simulations of somewha t more 
massive haloes (10 10 M©), | Governato et al. (|2010| found 



that increasing resolution was key in reducing the central an- 
gular momentum of galactic disks that form in their haloes 
and attribute this difference to better resolving spatially sep- 
arated clouds in their higher resolution simulations, which 
enhanced the effect of supernova feedback as a regulator of 
star formation. They also found convergence in behavior be- 
tween their medium resolution simulation (with a force res- 
olution of 116 pc) and their high resolution simulation (with 
a force resolution of 86 pc), in contrast to the behavior of our 
simulations R10 and R43. Therefore, while it appears that 
resolution in key in simulating star formation and super- 
nova feedback in galaxy simulations, the exact mechanisms 
at play and the scales on which behavior converge are de- 
pendent on the physical prescriptions and hydrodynamical 
methods used. 

Simulations R10 and R43 do not converge on a star- 
formation history or baryon fraction for our final halo, how- 
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ever we argue that many of the conclusions we have drawn 
from our higher resolution simulations hold. Clearly, reion- 
ization and feedback operate in the same way in both simu- 
lations, where the overall gas fraction is regulated by reion- 
ization and the dense gas evolution is regulated by supernova 
feedback. We do see an increasing effect of supernova feed- 
back on the gas fraction with increasing resolution, so it is 
possible that increasing the resolution even more will have 
an even greater effect. The effect of reionization, however, 
is unchanged and still quite dominant. We may not have 
converged on a final stellar mass for this halo. We therefore 
conclude that our measures of the final stellar mass is likely 
to be an upper limit. 



4.4 Metal ejection 

We find that the main mechanism for gas ejection from these 
low mass haloes is photo-evaporative winds due to reioniza- 
tion. It remains unclear, however, if this is also the mecha- 
nism for metal ejection. We see evidence for metal rich su- 
pernova heated winds in Figure [4] but it also appears from 
Figure [3] that the metallicity of dark haloes increases after 
reionization. 

We can answer this question by comparing simulations 
with and without supernova feedback and reionization: RIO, 
one of our canonical runs that contains our complete phys- 
ical model; RIO-noFB-LimCool, which has no thermal su- 
pernova feedback; and RIO-noUV, which has no UV back- 
ground. Figures [6] and [8] show that these three simulations 
have produced varying amounts of star particles over some- 
what different star formation histories, but we can neverthe- 
less draw robust conclusions. 

Figure [13] shows projections of gas metallicity of a large 
volume that encompasses the dwarf group at z = 3.3, which 
eventually forms the final dwarf halo. It clearly shows that 
in simulations with thermal supernova feedback, metals are 
ejected to much greater distances and more completely fill 
the volume. The simulation without thermal feedback, R10- 
noFB-LimCool, has significantly less IGM metal pollution. 
This is despite the fact that the total gas fractions of haloes 
in RIO-noFB-LimCool are quite similar to those seen in 
RIO (Fig. [8]). The gas fractions of haloes in RIO-noUV are 
quite elevated at this redshift as compared to haloes in RIO 
(Fig. [6]). Despite this, the extent of metal pollution is rel- 
atively similar between the two runs with supernova feed- 
back. Therefore, we conclude that supernova feedback is the 
primary mechanism for metal ejection in our model. 

As we have discussed, stellar populations of dSph galax- 
ies are extremely metal poor and this does not appear to be 
simply due to their low star formation rates. | Kir by, Mar-| 
tin, & Finlator ( 2011| have demonstrated that dSphs in the 
Milky Way eject almost all of their metals. Dwarf spheroidals 
are also observed to be quite gas poor ( |Mateo|1998| ). As we 
have discussed, environment most likely plays a key role in 
the gas content of Milky Way dwarfs. However, our simu- 
lations at least demonstrate that the main mechanism for 
gas expulsion need not be the same as for metal ejection. 
This picture is consistent with the simulations of Mac Low] 
fe Ferrara| ( |l999| ) and |Fragile et al.| p003) , who found, us- 
ing high-resolution, idealized models of dwarf galaxies, that 
metals were more easily ejected from galaxies than gas. 

The discrepancy in metallicity between our simulated 
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Figure 13. Gas density weighted projections of gas metallicity 
at redshift 3.31 for RIO, RIO-noFB and RIO-noUV. All panels 
show projections of 150 kpc cubic volumes that encompass all 
the progenitor haloes at this redshift. 
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stellar populations and observed systems described in sec- 
tion 4.1.2, indicate that our models are not pushing metals 
out effectively enough from star forming gas in haloes, so 
we cannot place any constraints on the density of metals in 
the IGM produced by dwarfs. It is certainly possible that 
feedback mechanisms other than supernovae may play an 
important role in ejecting metals, such as ionization from 
young stars. What our models do show, however, is that in 
systems where the gas fraction is determined by the cos- 
mic UV background, an extra ingredient that preferentially 
affects metal rich gas is needed to eject metals. 

If supernova feedback is primarily responsible for the 
stellar metallicity of dSphs, then these systems are an ex- 
cellent test to judge the efficacy of feedback models, given 
the clear constraints of observations. Several simulated mod- 
els produce low metallicities for Milky Way dSph ana- 
logues using a variety of approaches, such as idealized, non- 
cosmological models (e.g. |Revaz &; J ablonka 2012); cosmo- 
logical models focusing on high redshift evolution (|Ricotti 



|&~Gn edin 2005] |Tassis, Gnedin, fe Kravt sov 2012] |Wise et 



al.|2012bj ; coarsely resolved cosmological models of satellite 



systems of Milky Way type spiral galaxies (Okamoto et al. 
2010; Sawala, Scannapieco, & White 201 2]); and highly re - 



solved zoom-in simulations like our own ( |Sawala et al.|2010 ). 
Most of these studies rely on feedback schemes that compen- 
sate for resolutions too low to capture the Sedov phase of 



supernovae, with the exception of Wise et al. (2012b), who 



have very high resolutions (1 comoving pc) in their high 
redshift study and find that momentum from feedback, in 
the form of radiation pressure, is necessary to produce low 
metallicities. The model most comparable to our own, which 
is run down to z = 0, is that of I Sawala et al.l 120101) , who 



in high resolution simulations (< 100 pc) of isolated haloes 
similar in mass to our own, were able to produce haloes 
that lie on the observed luminosity-metallicity relation. It is 
unclear, however, how natural this result is given that the 
fraction of metals injected into the hot (or cold) phase can be 
calibrated in their subgrid model; in addition, we note that 
metal mixing remains a serious uncertainty in SPH ( |Booth| 
et al.|2012 ). Simulations of larger cosmological volumes that 
contain larger haloes have also found that feedback is im- 
portant in explaining the evolution of galaxy metallicity and 



enrichment of the IGM over cosmic time ( |Finlator &; Dav£ 
2008). 



4.5 Dark matter properties 

Dark matter haloes in cosmological N-body simulations fol- 
low a predictable density profile with a steeply-sloped cusp 
at the centre ( |Navarro, Frenk, &; White 1997 Navarro et 
al. |2010[) . Several previous studies have found that baryon 



physics, particularly supernova feedback, can weaken the 
central dark matter cusp (e.g. |Navarro, Eke, fe Frenk|[T996T 



Gnedin fc Zhao|2002| |Governato et al.|2012| . Several mech- 
anisms for this transformation have been proposed, such 
as the response of dark matter to the sudden removal of 



gas from halo centres (e.g. Navarro, Eke, & Frenk 1996 



Gnedin & Zhao 2002), and the direct transfer of energy 



from repeated energetic gas outflows to dark matter par- 



ticles ( [Pontzen &; Govern ato 20 12[ [Governato et al.||2012 ). 

Measuring the actual dark matter density profiles of 
galaxies is, of course, quite challenging, but one way is to use 




Figure 14. Dark matter density profile (solid black line), stellar 
mass density profile (dotted line) and gas density profile (dashed 
line) for the final halo at z = in R10. The dark matter density 
profile for the final halo in R10-DM is shown in grey. 



the resolved stellar populations in dwarf spheroidal galaxies 
To date, studies comparing the dynamics of dSphs to cosmo- 
logical N-body simulations have been inconclusive ( |Boylan- 
Kolchin, Bullock, & Kaplinghat 2012; Strigari, Frenk, &; 



White| |2010| |Vera-Ciro et al.||2012| ). Density profiles with 
shallow or steep central slopes have been found to fit the 



observational data equally well ( Walker et al.|2009b ). 

Figure [14] compares the dark matter density profiles for 
the final halo in R10 and R10-DM, a dark matter only sim- 
ulation run with the same initial conditions. We see that the 
dark matter density profile is almost the same between R10 
and R10-DM; the halo in R10 is slightly less dense in the 
inner regions relative to R10-DM, indicative of a mild heat- 
ing effect. We therefore conclude that the impact of baryon 
physics on the dark matter density was negligible for our 
halo. Although not shown in Figure |14| both profiles are 
well fit by a NFW profile. The stellar density profile does 
not follow a NFW profile. It appears to have a central slope 
similar to a NFW profile but falls off more quickly beyond 
300 pc. Visual inspection of the final stellar density map 
shown in the bottom right panel of Figure [4] also shows 
that dumpiness and tidal features from late mergers remain, 
which produces a bump in the stellar density profile. 

The one other high-resolution study to examine the den- 
sity profile of dark matter haloes at this halo mass scale with 
baryon physics is |Governato et aT] ( |2012| ), though using a dif- 
ferent hydrodynamical method and stellar feedback model. 
They find flattened central density profiles for haloes with 
virial masses above 10 10 M©. For haloes below this thresh- 
old, they find haloes down to our halo mass of 10 9 M have 
density profiles consistent with NFW profiles, similar to our 
results. Haloes in their models do not appear to have consis- 
tently flattened central density profiles until the stellar mass 
rises above 5 x 10 6 M . 

As we have discussed, our supernova driven outflows 
may not be sufficiently energetic (as evidenced by the metal- 
licity distribution). It is possible that this failure results in 
the central density of our halo remaining cuspy, however it 



© 2012 RAS, MNRAS 000, [l]-?? 



The effect of feedback and reionization on star formation in low-mass dwarf galaxy haloes 23 



is also possible that the low star formation rates of our halo 
are just insufficient to produce enough energy to have a sig- 
nificant impact on the density profile of our simulated halo. 



4.6 Neglected physics 

In this section, we briefly discuss physical effects that we 
do not include in our model. We do not have an implemen- 
tation of local UV heating from young stars, either ioniz- 
ing or dissociating, in our model. This effect may have an 
impact on dense gas as young stars will be located within 



close proximity. |Gnedin ( 2010| found that local UV heat- 
ing sources dominated over the cosmic UV background for 
both ionizing and Lyman- Werner radiation in haloes with 
masses between 10 10 and 10 12 Mq at a redshift of three. We 
are probing haloes at much lower masses with much lower 
star formation rates, but it may be that UV heating from 
young stars can suppress star formation or make supernova 
feedback more effective by preheating gas. 

We do not include any effects of dust in our model. The 
formation of H2 on the surface of dust grains is impo rtant in 



gas with metallicities above approximately 10 Z (Glover 



2003). The main effect of H2 in our models is as a coolant in 



very metal poor gas. Since the amount of dust in gas with 
metallicities below 10 -3 Zq is too small to be important in 
H2 production, the lack of dust does not significantly affect 
cooling and therefore star formation in our models. It does 
mean that our models do not accurately follow the molecular 
gas fraction of haloes, which may be an interesting topic of 
future study. 

We do not model early metal-free star formation as a 
distinct mode of star formation. |Wise et al.| ( |2012a| pre- 
sented high-resolution simulations of low mass haloes at high 
redshift that include the transition between Population III 
and Population II star formation. They found that a single 
Population III supernova can enrich the ISM of a halo to 
the level where metal line cooling dominates over molecular 
H2 cooling. 

We model supernovae as producing exclusively ther- 
mal energy. Supernovae also produce substantial amounts 
of radiation, which can not only heat the gas, but can also 
exert radiation pressure. The computational cost necessary 
to do radiative transfer in cosmological calculations is cur- 
rently too expensive to follow the evolution of haloes down 
to z = 0, however several groups have explored its effect at 
high redshift (e.g. Wise et al. 2012b). There are also other 
effects such as magnetic fields and cosmic rays which we ne- 
glect and may potentially be important in galaxies of this 
size ( jWadepuhl fe Springel||20TT| . 



5 CONCLUSIONS 

We have conducted the highest resolution cosmological sim- 
ulations to date of a 10 9 Mq dwarf halo down to a redshift 
of zero. Our resolution of 11 comoving pc means we achieve 
physical resolutions of a few parsecs during the epoch of 
greatest star formation in our haloes, between redshifts of 6 
and 10. Our canonical runs includes metal cooling, molec- 
ular hydrogen formation and cooling, photoionization and 
photodissociation from a metagalactic background (with a 
simple prescription for self-shielding), star formation, and a 



model for supernovae driven energetic feedback. Our main 
results are summarized below: 

(i) We find that reionization is primarily responsible for 
setting the total gas fraction of our haloes. This occurs 
through photo-evaporative winds, which expel most of the 
gas envelope in our simulated systems (although not the cen- 
tral, dense core). Once the gas is ejected, our halo is unable 
to efficiently re- accrete new baryons, and so the gas fraction 
remains suppressed to z = 0. 

(ii) The dense gas, from which star particles form, ap- 
pears to be mostly regulated by feedback from star forma- 
tion, although UV heating contributes. 

(hi) The timing of reionization has a major impact on the 
final stellar mass of our halo, by interrupting the cooling of 
gas onto lower mass progenitor haloes. Depending on the 
patchiness of reionization, an early or delayed epoch of ion- 
ization may result in an order of magnitude change to the 
resulting stellar mass 

(iv) Our final halo has a mass-to-light ratio consistent 
with Milky Way dSphs. The systems our models most re- 
semble are the ultrafaint dSphs found in SDSS, although, 
our simulated halo is located in a much lower density envi- 
ronment than the observed ultrafaints. The star particles in 
our models are generally too metal rich compared to Milky 
Way dSphs with similar luminosities, indicating a problem 
with our feedback model's ability to eject metals. Finally, 
our isolated halo is extremely gas poor at z = 0, in contrast 
to many dwarfs in the Local Group and the local volume 
in low density environments but in agreement with many of 
the ultrafaint dSphs. 

(v) While reionization is the main mechanism responsible 
for expelling gas from progenitor haloes, supernova feedback 
is mostly responsible for expelling metals. 

(vi) H2 line cooling is a crucial and necessary ingredient 
to start star formation in these low mass haloes. Atomic line 
cooling in the absence of metals is insufficient by itself to 
cool gas to the densities required for star formation in such 
shallow potential wells. Without a first generation of star 
formation to pollute the ISM with metals, there is no coolant 
sufficiently effective to cool gas at subsequent epochs. 

(vii) We find that the maximum star formation rate and 
overall duration of star formation is sensitive to the spatial 
resolution of the calculation. Our models appear to have 
not yet converged to a star formation history for our halo, 
but we argue that the final stellar masses we find for our 
simulated halo are upper limits. 

(viii) Our final halo, which has a mass of 10 9 Mq at z = 0, 
forms hierarchically from several progenitor haloes, many of 
which are star-forming. The stellar population in the final 
halo was therefore built up from star formation episodes in 
these separate haloes. 

We conclude that while we do reproduce many of the 
properties of dwarf systems, such as high mass-to-light ra- 
tios and low luminosities, more work is required to under- 
stand how stellar feedback works to regulate star formation 
and metal ejection in low mass dwarf galaxies, and hence 
the properties of stellar populations in dwarf haloes. The 
feedback scheme we adopted was deliberately simple, and 
not tuned to reproduce observations. Future work will in- 
volve more realistic modeling of supernovae, a model of local 
sources of ionization, an expansion of the halo we simulate 
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to larger masses and potentially haloes in different environ- 
ments. 
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